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A short summary of results of recent simulations of (un)polarized Drell-Yan exper- 
iments is presented here. Dilepton production in pp, pp, 7r~p and 7r+p scattering 
is considered, for several kinematics corresponding to interesting regions for exper- 
iments at GSI, CERN-Compass and RHIC. A table of integrated cross sections, 
and a set of estimated error bars on measurements of azimuthal asymmetries (as- 
sociated with collection of 5, 20 or 80 Kevents) are reported. 



1. Introduction 

The aim of this work is to give some useful reference numbers for 
planning Drell-Yan experiments aimed at the measurement of transverse 
spin/momentum related azimuthal asymmetries. 



Table 1. Total cross sections (nb) for several colliding particle combinations, mass 
ranges and S- values. 



M (GeV/c2) 


1.5-2.5 


4-9 


12-40 


1.5-2.5 


4-9 


12-40 


S (GcV2) 




PP 






PP 




30 


0.03 


<<pb 


<<pb 


1.3 


0.3 pb 


<<pb 


200 


0.7 


0.01 


<<pb 


4.4 


0.35 


<<pb 


(200)2 


3.5-17 


1.2-5.5 


0.06-0.26 


5-18 


2.6-7 


0.24-0.4 












7r+p 





30 0.9 1 pb <<pb 0.25 0.1 pb <<pb 

200 1.9 0.25 <<pb 0.7 0.07 <<pb 

(200)2 1.8-5.6 0.75-2.1 0.1-0.2 1.5-4.8 0.5-1.4 0.04-0.1 

Note: For the high-energy case, two different parameterizations have been used, 
leading to pairs of a values. See text for details. 
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Figure 1. pp event scatter plot for 5 = 200 GeV^, and dilepton mass in the range 4-9 
GeV/c2. 



This includes total cross sections for several kinematical options for pp, 
pp and n±p Drell-Yan dilepton production: squared hadron-CM energy S 
= 30 GeV^ 200 GeV^, (200)^ GeV^ and dilepton masses in the ranges 
1.5-2.5 GeV/c2, 4-9 GcV/c^, 12-40 GeV/c^. For some relevant situations, 
estimates of the asymmetry error bars are reported, for sets of 5, 20, 80 
Kevents divided into 10 bins of the longitudinal fraction x. Previous Drell- 
Yan data and fitting relations^^^ are the basis of the initial core of the used 
simulation code. Several details on the formalism, together with the most 
recent examples of simulated asymmetries, are presented elsewhere in this 
workshop^, and in published work by myself and M.Radici^'^. 



2. Total Cross Sections. 

Total cross sections are shown in table 1. For the two lower S'— values 
they have been evaluated with the difl:erential cross-section fit relations^ 
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Figure 2. pp event scatter plot for 5 = 30 GeV^, and dilepton mass in the range 1.5-2.5 
GeV/c2. 



coming from measurements of 7r~A and pA at S = 250-400 GeV^. The 

two cases Tr~p and pp correspond to these cross sections for Z/A = 1. 
For TT+p we have assumed n~p = n^p for the pion sea contribution, 
'K~p = (l/4)7r"'"n for the pion valence contribution, and calculated 'jt~A 
for Z/A = 0. Cross sections for pp have been evaluated by substitut- 
ing the {sea + valence){sea + valence) structure of pp with the struc- 
ture l/2[sea(sea + valence) + {sea + valence)sea]. For the largest S case 
the calculation based on the previous parameterization (with sea ^ const 
for small x) has been sided by an alternative calculation using the (LO- 
intermediate gluon) MRST distributions ^, with sea ~ x~'^ (A « 0.2-0.3 for 
mass < 9 GeV/c), and based on data sets including several recent Drell- 
Yan meaasurements^. No evolution was applied, which is unproper for M 
» 10 GcV/c^. The ii'— factors have been assumed as constant but tuned 
to reproduce with both methods the measured cross sections at S = 250 
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Figure 3. pp event scatter plot for S = (200)^ GeV^, and dilepton mass in the double 
range 4-9 and 12-40 GeV/c^. 



GeV^, masses 4-9 GeV/c^^. The former strategy leads to the smaller re- 
ported cross section values, the latter to the bigger ones. For the w^p 
high-energy case double values refer to different parameterizations for the 
distribution functions of the proton only, pion distribution functions have 
not been changed. For the lowest mass range 1.5-2.5 the smaller number 
of contributing quarks introduces a reduction factor w 1/2. The difference 
between the other two mass ranges is not essential. 



3. Event distributions 

The event distribution has the general form 



NiS,x,x',PT,0 = F{S,x,x',Pt)-[1 + A{x,x',Pt,0]- W 
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Figure 4. X2- integrated distribution of the scatter plot of fig.l. Events are divided into 
50 Xl-bins. 



where x,x' , Pt describe the virtual photon kinematics in the hadron center 
of mass {Pl/{S/2) = x — x' , M'^ / S = xx'), while ^ represents compactly 
the set of variables describing the angular distribution of the leptons in the 
CoUins-Soper frame (= the photon polarization). F{S, x, x\ Pt) alone gives 
the virtual photon event distribution. A{x , x' , Pt , £,) averages to zero over 
all the solid angle, and describes the asymmetry properties of the lepton 
distribution in unpolarizcd, single or double polarized DY. 

The scatter plots of figs 1,2 and 3 reproduce the event distribution 

N{S,x,x') = J (fPTFiS,x,x',PT) (2) 

for some relevant pp and pp cases. Fig. 4 reports 

N{S,x) = J dx'N{S,x,x') (3) 

where the integrated distribution is the one of fig.l. Events in the xx' 
scatter plots concentrate near the hyperbole xx' = M'^rnin/S, because a 
(X q{x)q{x')/A'P . For the same reason, in the case of fig. 3 the lower event 
band (M in the range 4-9 GeV/c^) contains 95 % of all the events reported 
in the figure. 
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asymmetry error bars for 5, 20, 80 Kevents 
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Figure 5. Error bars on the azimuthal asymmetry for pp, S = 30 GeV^, dilepton mass 
in the range 1.5-2.5 GeV/c^. 



4. Asymmetry Error Bars 

The integrated N{S,x) distribution has its peak at a; « 3M'^min/S. To 
the right of the peak, N{S, x) ~ l/x" with n > 1. This selects for each 5, 
Mmin the a;— range where error bars are smaller. The shown error bars in 
figs. 5, 6 and 7 refer to the special case of Sivers asymmetry, however they 
are the same for any kind of left/right asymmetry with respect to a CoUins- 
Soper azimuthal angle 4>. The asymmetry is defined as {A — B)/{A + B), 
where A and B arc the event numbers with positive or negative sin{(h — (j)s). 
Error bars have been calculated by assuming constant 0.05 asymmetry ev- 
erywhere, and repeating the simulation 10 times to calculate fluctuations. 
Despite the size of the error bars is reasonably stable for N repetitions > 5, 
with N repetitions = 10 there are still small but evident fluctuations in the 
error bar size. This is due to most x— bins being fllled with event numbers 
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Figure 6. Error bars on the azimuthal asymmetry for n p, pion beam energy 100 GeV, 
fixed target, dilepton mass in the range 4-9 GeV/c^. 

N{+) + N{-) < 1000. For iV(+) + N{-) < 50 error bars arc not shown. 
Typically the largest shown error bars refer to event numbers N{+) + N{—) 
~ 100. Examination of the error bar fluctuations suggests that the reported 
error bars are reliable within a factor \f2. If one repeats the calcu- 

lation by assuming asymmetry or 0.1, systematic changes of these (purely 
statistic) error bars are smaller than the above fluctuations. So, unless they 
damp asymmetries by orders, all those reducing coefficient like polarization 
dilution etc are not influent on the error bar sizes. For asymmetry > 33 
%, A'^(— ) < N{-\-)l2. So, in the case of really large asymmetries, the con- 
sequent small value of N[—) in those bins whose population is << 1000 
introduces errors whose size may be uncontrollably larger than the esti- 
mated ones. In this case the error simulation must be more specific to be 
appropriate. 
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asymmetry error bars for 5, 20, 80 Kevents 
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Figure 7. Error bars on the azimuthal asymmetry for pp, S = (200)^ GeV^, dilepton 
mass in the range 12-40 GeV/c^. 
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